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Abstract 

We report new measurements of the aeoustie excitation of an A15056 su- 



<N 

perconductive bar when hit by an electron beam, in a previously unexplored 

temperature range, down to 0.35 K. These data, analyzed together with previ- 
ew ous results of the RAP experiment obtaincxl for T > 0.54 K, show a vibrational 

response enhanced by a factor ~ 4.9 with respect to that measured in the nor- 
l/^ mal state. This enhancement explains the anomalous large signals due to cosmic 

rays previously detected in the NAUTILUS gravitational wave detector. 

■T-H Keywords: Gravitational wave detectors, Cosmic rays. Radiation acoustics 

^ PACS: 04.80.Nn, 41.75.Ht, 96.50.S-, 61.82.Bg, 74.70.Ad 

• ^ 

^ 1. Introduction 

Cosmic ray showers can excite sudden mechanical vibrations in a metallic 
cylinder at its resonance frequencies; in experiments searching for gravitational 
waves {gw) these disturbances are hardly distinguishable from the searched 
signal and represent an undesired source of accidental events, thus increasing 
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the background. This effect was suggested many years ago and a first search was 
carried out with Hmited sensitivity with room temperature Weber type resonant 
bar detectors and ended with a null result [1] . Later on, the cryogenic resonant 
gw detector NAUTILUS [5] was equipped with a streamer tube extensive air 
shower detector [3] and the interaction of cosmic ray with the antenna has been 
studied in detail. This apparatus allowed the first detection of cosmic ray signals 
in a gw antenna, that took place in 1998, when NAUTILUS was operating at a 
temperature T = 0.14 K [4], i.e. below the superconducting (s) transition critical 
temperature Tc ~ 0.9 K. During this run many events of very large amplitude 
were detected. This unexpected result prompted the construction, in 2002, of 
a scintillator cosmic ray detector also for the EXPLORER gw detector as well 
as the beginning of a dedicated experiment (RAP) [S], that was planned at the 
INFN Frascati National Laboratory to study the vibration amplitude of a small 
A15056 bar caused by the hits of a 510 MeV electron beam. The experiment 
was also motivated by the need of a better definition of the thermophysical 
parameters of the alloy A15056, used in the bar detector, at low temperatures. 
A detailed study of this effect is indeed useful to study the performance of 
gw bar detectors for exotic particles [B] and to understand the noise due to 
cosmic rays in interferometric gw detectors [7]. In this paper we summarize 
our previous knowledge on this effect. We then report the final results of the 
RAP experiments presenting measurements down to 0.35 K, and show how 
these new data help in sheding light on the 1998 anomalous NAUTILUS high 
energy events. We also recall that the detailed study of this thermo-acoustic 
effect has applications in devices used to monitor and measure particle beams 
characteristic, and in particular in monitoring high power beams |8]. 

2. The Thermo-Acoustic Model 

2.1. The model: normal conductive (n) state 

B.L. Beron and R. Hofstadter OIlO] first measured mechanical oscillations in 
piezoelectric disks hit by a high energy electron beam. The authors first pointed 
out that cosmic ray events could excite mechanical vibrations in a gw metallic 
antenna and that, consequently, cosmic rays could represent a background for 
experiments aimed at the detection of gw. The interaction of a ionizing particle 
with the bulk of a suspended cylindrical bar generates a pressure pulse in the 
bar. More in detail, the energy lost by the particle in the bar causes a local 
warming up of the material; the local thermal expansion in the bulk generates 
the pressure wave. This sonic pulse determines the excitation of the vibrational 
elastic modes of the suspended bar. 

A.M. Grassi Strini et al. |TT] reported the results of an experiment based on 
a pure aluminum bar exposed to a proton beam. The experimental data were 
compared to a theoretical model based on the Fourier response of a thin bar to 
the pressure wave originated by a delta-like thermal perturbation. If the heating 
is in the bar center, only even Fourier harmonics are allowed. The "maximum 
amplitude of oscillation for the fundamental longitudinal elastic mode" , in the 
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following referred to as "Amplitude", for a material in a normal (n) state of 
conducting is given, according to this Thermo- Acoustic Model (TAM), by: 

where the suffix "th" stand for the theoretically expected value. This result 
applies to a thin cylinder (with radius R and length L, R <^ L and mass M), 
for a beam hitting on center of the cylinder lateral surface. Here W is the 
total energy released by the beam to the bar, a is the linear thermal expansion 
coefficient and cy is the isochoric specific heat. The dimensionless Griineisen 
parameter 7 of the material includes the a/cy ratio: 

(2) 

pcv 

where /3 is the volume thermal expansion coefficient (/? = Sa for cubic elements), 
Kt is the isothermal bulk elastic modulus and p is the material density. The 
Griineisen 7 slightly depends on the temperature when the material is in the n 
state. 

Eqn. ([T]) is a limit case of a more general problem, when the paths of the 
interacting particles in the bulk [121 lEl HI] are considered. Introducing a vector 
field u(x, t) describing the local displacements from equilibrium, the amplitude 
of the k — th mode of the cylinder oscillation is proportional to: 
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(3) 



where AP*'"^''™ is the pressure pulse due to the sonic source described above, 
dW/dx is the specific energy loss of the interacting particle. A' is the cross 
section of the tubular zone centered on the particle path in which the effects 
are generated and Tk = J dl{\7 ■ Ufe(x)) is a line integral over the particle path 
involving the normal mode of oscillation u/j(x). The Amplitude predicted by 
Eqn. ([1]), can be rederived from Eqn. Q in the simplified case of a thin bar 
{R/L ^ 1) and for a particle hitting on the bar center. We can therefore adjust 
the value of Amplitude predicted by Eqn. ([!]) to a more correct value: 

where e is a corrective parameter estimated by a Monte Carlo (MC) simula- 
tion [S], which takes into account the solutions 0[(i?/i)^] for the modes of 
oscillation of a cylinder, the transverse dimension of the beam at the impact 
point and the trajectories of the secondary particles generated in the bar. The 
value of e for the bar used in the experiment is estimated by MC to be -0.04. 
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2.2. The model: superconducting (s) state 

When the material is in the s state, a part of the energy lost by particle 
causes the suppression of the superconductivity in a region, called hot spot, 
that is centered around the particle path. The maximum possible radius of the 
hot spot, vhs: is obtained by equating the specific energy lost by the particle, 
dW/dx, to the enthalpy variation (per unit volume), Ah, for the transition from 
the s state at temperature T to the n state [HI [161 113 HH] : 




TT Ah 



(5) 



where A" is the cross section of the zone switched to the n state. 
The creation of a hot spot by a particle interacting with a material in s state 
causes a further correction to Amplitude of Eqn. Q, a term which is peculiar to 
the particle propagating in a zone now switched to the n state. The additional 
contribution to the amplitude of the cylinder oscillation mode k is proportional 
to [ni[T3|: 
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Here AV^ and AS are the differences of the volume and entropy in the two states 
of conduction. The differences can be expressed in terms of the thermodynamic 
critical field He and it follows [TU], in first approximation, thai[^ 

AV ^ Vn-Vs ^ HcdHc 
V V 4tt dP 

AS ^ Sn- Ss ^ dH, 
V V A-K dT 

The quadratic dependence Hc{t) — Hc{0){l — t^) on t, where t = T/Tc, is 
assumed in computing the differences. 

Therefore the value of Amplitude due to a particle creating hot spots in a 
material in s state is given by: 
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and the definition ([2]) of 7 is used to obtain: 



(6) 



•^We keep the CGS electro-magnetic system of units for the magnetic field, as used by the 
authors of the cited articles, and we convert the density of the magnetic energy to SI units. 
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T 


a 


cy 




[K] 


[10-6 K-i] 


[J mol-i K-i] 


[10-1*' m J-i] 


264 


22.2 


23.5 


2.23 


71 


7.5 


7.94 


2.23 


4.5 


5.8 X 10-3 


7.6 X 10-3 


1.80 


1.5 


1.5 X 10-3 


2.1 X 10-3t 


1.72 



Table 1: Normal state of conduction. Amplitude normalized to the beam released energy 
W and input values for the calculation (a,CY) in the case of the RAP bar (L = 0.5 m; 
M = 34.1 kg) made of pure aluminum. The correction of Eqn. Q is applied, ^cy value for 
A15056. 
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Eqn. Q and ^ show, by inspection, that the Amphtude B hnearly de- 
pends on W, the energy released by particle, both in the n and in the s states. 
Therefore it appears natural to consider, as we do in the following, the ratio 
B *''/W as a measure of the relevant material properties. Finally we note that 
the knowledge of the specific heat of the material for the s state, c^, allows us to 
approximate the exact Eqn. ([5| with the following condition for the transition 
s — >■ n of a volume V of the material at temperature T, due to the absorption 
of energy W from the particle [151 HZ] : 

W > VCi{T) (7) 

with: ^ 

Ci{T) = I \s{T')dT' 
Jt 

Moreover, the knowledge of Cs{T) allows to derive Ah from the relation Ah{T) = 
C/(T)+T^ [20]. 

2.3. Amplitude predictions for the normal and superconductive state 

In order to compute, by means of Eqn. Q, the expected value of Amplitude 
at different temperatures, we need, for the n state, both a{T) and cy(T) of 
the material. As these values are not well known for the A15056 alloy, we 
used those of pure aluminum. Polynomial interpolations on data of Ref. [21] 
(12 < T < 300 K) and the parametrization in Ref. ^ (T < 12 K) give a(T), 
while cv{T) is obtained by polynomial interpolations on values of cp reported 
in Ref. [53]. Table [T] shows the computed values of a, cy and the normalized 
Amplitude B„*''/W. 

For the s state instead, measurements performed at very low temperatures 
[24| on samples belonging to the same production batch of our A15056 bars, allow 
us to characterize the relevant properties of that alloy. The measurement of the 
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transition temperature to the s state, carried out using the mutual inductance 
method, yields Tc = 0.845± 0.002 K and a total transition width of about 0.1 K. 
Few data of specific heat for A15056 are available in the literature [25] ; however, 
new cy measurements were performed above and below |24) . 

In the temperature interval 0.9 < T < 1.5 K, i.e. in the n state, assuming for 
the specific heat the usual low temperature parametrization: cy = TT + "^T^, 
the measurements give the values F = 1157 ± 31 erg cm^'^ K^^ for the elec- 
tronic specific heat coefficient and 4' ~ 140 ± 10 erg cm"'' K"'' for the lattice 
contribution. 

On the other hand, computing the Amplitude in the s state by means of 
Eqn. ^ requires the knowledge of a) the thermophysical parameters a„ and 
cv,n of the material, in order to evaluate Bn for the n state below Tc and b) 
the dependence of on T and P for calculating the derivatives OHc/dT and 
dHc/dP. Fhe requirement a) cannot be fulfilled due to the lack of knowledge 
of a„ for A15056 and we therefore assume for B^/W the value measured just 
above T^- This assumption is justified by the fact that 7„ usually has, below Tc, 
a very weak dependence on temperature. Regarding requirement b), we derive 
dHc/dT a,t T < Tc from the He parabolic dependence on t; we also assume that 
the unknown dependence of dHc/dP on t at P = for A15056 is equal to that of 
pure aluminum and, therefore, can be obtained from the tabulation of He as a 
function of T and P contained in Ref . . If the superconducting properties of 
A15056 can be described by the BCS theory, then iJc(O) « 2.42 F^/^t; ~ 70 Oe. 
Insertion of the numerical values in the Eqn. ^ yields to values of Bg/W 
ranging from -9.2 x 10"^*^ to -7.3 x 10"^'' m J in the temperature interval 
having limits 0.3 and 0.8 K, respectively. The lower limit of the temperature 
interval is constrained by the data availability in the Hc{P,T) tabulation of 
Ref. 

3. The RAP experimental setup 

3.1. The bar and the piezoelectric ceramics 

The RAP experiment has been fully described in Ref. fSl . Here we briefly re- 
call that the test mass is a cylindrical bar {R = 0.091 m, L = 0.5 m, M ~ 34.1 kg) 
made of A15056, the same aluminum alloy (nominal composition 5.2% Mg and 
0.1% of both Cr and Mn) used for NAUTILUS. The bar hangs from the cryostat 
top by means of a multi-stage suspension system ensuring attenuation from the 
external mechanical noise of -150 dB in the 1700-6500 Hz frequency window. 
The frequency of the fundamental longitudinal mode of oscillation of the bar is 
/o = 5414.31 Hz below T = 4 K. 

At temperatures below 10 K, the A15056 intrinsic Q factor is 4.1 x 10'' [27] 
corresponding to a decay time of the order of 20 minutes: a shorter decay time 
is desirable in order to have a more manageable repetition rate of the hits. This 
is simply achieved by the presence of a thermometer on one of the bar end 
face, that damps the oscillations to a decay time of the order of 30 seconds. 
Therefore, we only have to wait a couple of minutes between consecutive hits 
to avoid the pile-up of the signals. 
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Two piezoelectric ceramics (Pz), electrically connected in parallel, are in- 
serted in a slot milled out in the center section of the bar, opposite to the bar 
suspension point, and are squeezed when the bar shrinks. In this Pz arrange- 
ment the strain measured at the bar center is proportional to the displacement 
of the bar end faces. The Pz output is first amplified, and then sampled at 
100 kHz by an ADC embedded in a VME system, hosting the data acquisition 
system. A band pass filter between 300 Hz and 50 kHz is used to reduce the 
low frequency power line noise and to avoid Fourier aliasing. 

The measurement of the Pz conversion factor A, relating voltage to oscillation 
amplitude, is accomplished according to a procedure [55] based on the injection 
in the Pz of a sinusoidal waveform of known amplitude, with frequency /o and 
time duration less than the decay time of the mechanical excitations and on the 
subsequent measurement of Amplitude. The procedure is correct ii R/L ^ 1 
and carries a 6% systematic error. 

The value of A during the 2009 run was 1.26 x 10^ V m^^ at room temper- 
ature and 1.16 x 10^ V m~^ constant at temperature T < 4.5 K. 

3.2. The Frascati DA<^NE Beam Test Facility 

The DA'I'NE BTF transfer line can transport and deliver in controlled way 
electron or positron beams, from the end of the DA$NE Linac to a 100 to? 
experimental hall, where users normally set up and operate their tests and 
experiments. Particles can be provided in 1 or 10 ns duration pulses, with an 
injection frequencjj^ spanning from 1 to 49 Hz, in a wide energy range (25 - 
750 MeV for electrons) and intensity (from 1 up to 10^° particles/pulse). 

Since the end of 2004, when the RAP experiment was first operated with 
the aluminum bar, some important upgrades [29' and diagnostic improvements 
|30j have been accomplished. The most important improvement concerned the 
installation, during the 2006 shutdown, of a pulsed dipole magnet at the end 
of the Linac: such magnet allows to alternate the beam between the DA<I>NE 
damping rings and the test beam area, thus enhancing the BTF duty cycle, 
due to the reduced switching time. So that, when beams are injected into 
the collider, the BTF can still receive beam, although with a lower repetition 
rate. The continuous feeding of the BTF line, during the time sharing with the 
DA$NE collider, is technically possible because not all the Linac bunches are 
needed for filling the accumulator rings. Obviously, in this operation scheme, 
the pulse length and the primary beam energy are the same as those fed to the 
DA$NE collider. Nevertheless this does not constitute a real limitation, since 
the facility is mainly operated in single particle mode (electrons/positrons), 
which is the ideal configuration for detector calibration and testing. In this 
case the beam characteristics are mainly determined by the magnetic field of 
the energy selector (the upstream 45° dipole, on the BTF transfer line) and by 
scrapers suitably positioned along the line itself. 



^Actually one of the 50 pulses in 1 s is sent to an hodoscope in order to reconstruct the 
energy profile of the beam at the end of the Linac 
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Figure 1: Schematic layout of the RAP setup in the BTF hall. 



The RAP experimental set-up was installed about 2.5 m away from the exit 
of the beam (see Fig. [T]) . The correct position of the beam at the exit of the 
line and at the entrance of the cryostat was monitored shot-by-shot by two 
high sensitivity fluorescence flags, 1 mm thick alumina doped with chromium 
targets. The spot size at the entrance of the cryostat, 50 cm far from the 
center of the bar, is about 2 cm in diameter. The Monte Carlo simulation of 
the detector indicates that the beam spot at the surface of the cylindrical bar 
preserves almost unchanged, due to the negligible effects of the cryostat vacuum 
and thermal shields intercepted by the beam. 

The beam multiplicity for the RAP measurements performed in 2009 spanned 
from 10^ to 10^ electrons per bunch of 10 ns width. The current intensity mod- 
ulation was obtained by properly changing the aperture of the tungsten slits 
(both horizontal and vertical), along the BTF transfer line (SLTBij in Fig. [2]): 
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Figure 2: Magnetic layout of the BTF transfer line: DHSTPOOl is the pulsed dipole magnet 
added during 2006 upgrade. 



this procedure, unlike that used in the previous RAP measurement campaign 
and based on the defocusing of the beam by quadrupoles, allows us to avoid the 
beam degradation, since the particles in the external tails are cut away. 

For the beam diagnostics, a monitor (WCM), based on an integrating cur- 
rent transformer, readout by a charge digitizer, was used to measure the pulse 
charge. The device accuracy is 3% and the readout noise fluctuations give a 
measurement error a — 1.5 x 10^ electrons. The device is equipped with a 
calibration coil used to control the gain, the noise and the time shaping of the 
generated signals. 

3.3. Cryogenic setup 

The RAP cryogenic setup consists of a KADEL commercial liquid helium 
cryostat, 3.2 m high and 1 ni in diameter, suspended on a vertically movable 
structure, and containing a dilution refrigerator. A schematic view of the cryo- 
stat together with the cold side of the dilution refrigerator is depicted in Fig. [3j 
The liquid helium (LHe) and liquid nitrogen (LN2) dewars, with a capacity of 
340 L and 200 L respectively, are placed in the upper half. Three stainless steel 
cables are suspended from the top flange to support the experimental appara- 
tus. To avoid the radiation input, 8 aluminum radiation shields are mounted 
between the top flange at room temperature, the 77 K OFHC (Oxygen Free, 
High Conductivity) copper flange and the 4.2 K OFHC copper flange. These 
two latter flanges are mechanically connected with the LN2 and LHe dewars, 
respectively. The experimental chamber is positioned on the lower half of the 
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cryostat and is surrounded by one OFHC copper radiation shield connected to 
the Still flange, two aluminum radiation shields connected to the LHe dewar 
and the LN2 dewar respectively, and the outer aluminum container. The LHe 
container is indium sealed to separate the experimental chamber volume from 
the insulation zone. 

The cryostat hosts a continuous flow, closed cycle ^He-''He dilution refriger- 
ator, made by Leiden Cryogenics, with a base temperature of about 100 mK and 
a cooling power of about 1 mW at 120 mK. The continuous flow is ensured by a 
pumping system, composed of two Varian TV 551 NAV turbo-molecular pumps 
and an Edwards XDS 35 scroll pump. A Gas Handling System control panel 
manages the mixture flow in the circuit lines, either automatically or manually. 
In automatic operation a CPU running a software program, that reads mea- 
sures from several Pirani pressure gauges and a flow meter placed in the circuit 
line, manages the flow operating a number of solenoid valves. Temperatures 
of the experimental setup are measured by 11 thermometers of three different 
types (PtlOOO, FeRh and Ru02 resistances), connected to an AVS-45 and an 
AVS-47 Picowatt Resistance Bridges. The resolution of the RUO2 thermome- 
ters at the lowest temperatures is 0.25 mK. The cryostat is also equipped by 
4 vacuum gauges and a LHe level gauge. Inside the dilution refrigerator there 
are 2 capacitance gauges that measure the liquid level in the IK Pot and in the 
Still. All diagnostic data are gathered, via serial and GPIB interfaces, by a PC 
running a Lab View program which displays the readings on a synoptic window 
and records all the measurements. 

To avoid transmission of mechanical vibrations to the bar, the thermal links 
between the cold spots of the refrigerator and the experiment are kept to a 
minimum (n. 8 in Fig. [3|: i) a couple of soft, thin OFHC copper sheets between 
the Mixing Chamber and the top of the suspension (n. 9 in Fig. [s]), and ii) 
3 sheets, same as above, between the Still flange and the 0.6 K flange (n. 6 
in Fig. [3]). These contacts assure the bar and suspension cooling when the 
dilution refrigerator is in operation, i.e. below 4 K. Above this temperature 
gas conduction provides the heat removal by inserting a few mbar of gaseous 
helium in the experimental chamber. The gas is then removed before reaching 
the liquefaction (about 5 K). 

Cooling the cryostat to 80 K and filling the LN2 dewar takes about 3 days 
and about 1000 L of LN2. 12 hours and about 800 L of LHe are sufficient to 
cool the system from 80 K to 4.2 K and leave about 150 L of liquid in the dewar. 
The cryostat consumption, once thermalized, is about 1 L/h of LN2 and about 
1.5 L/h of LHe, that raises to about 2 L/h when the IK Pot is in operation. 

The minimum temperature reached by the bar during the data taking was 
343 mK, read by a RUO2 thermometer placed at the center of one of the bar 
end faces. 

3.4- RAP data collection and analysis 

The piezoelectric signals were recorded by an ADC sampled at 100 kHz run- 
ning under Lab View control. At low temperature, due to the low specific heat, 
the thermometer has sufficient sensitivity to measure the increase of the bar 
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Figure 3: Schematic view of the RAP cryostat: [1] LN2 reservoir, [2] LHe container, [3] 
Stainless Steel (SS) suspension cables, [4] 77 K flange, [5] 4-^ K flange, [6] 0.6 K flange, 
[7] SS screw with Teflon ring, [8] soft; copper thermal contacts, [9] copper suspension, [10] 
bar, [11] dilution refrigerator cold end, [12] IK Pot, [13] radiation shields. The filled squares 
represent the thermometers. 
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Figure 4: Bar temperature measured during a run. Four beam shots are visible. The steps in 
temperature are proportional to the energy deposited in the bar. They can be used to measure 
the A15056 specific heat, yielding a result in agreement with the measurement in Ref. 
Moreover, this information is useful to monitor the beam and to check the intensity measured 
by the current monitor. 



temperature after each shot. This additional information helped in monitoring 
the beam intensity. It was important to read this raise immediately after the 
hit, as the temperature rapidly relaxes back to its equilibrium value. For this 
reason, after the cooldown and during the beam measurement, we gave up mul- 
tiplexing the thermometers, so that the thermometer on the bar end face can 
be read out with a rate of 1 Hz. 

Fig. [4] shows the bar temperature measured during a low temperature run: 
it can be clearly seen that the bar is warmed up by the beam, and after each 
hit it relaxes to a higher temperature. To deal with this problem, the measure- 
ments were started at base temperature (343 mK) and taken with successively 
increasing temperature. 

The 100 kHz ADC data were processed both online for a fast response and 
offline with a more sophisticated procedure. Fig. [5] and Fig. |6] show several 
important features of the signal. Fig. [5] shows the ADC output for a typical 
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Figure 5: Above: an example of a shot at T = 0.368 K. After the shot the temperature raised 
to 0.405 K. The deposited energy was 30 mj. The horizontal units are samples, the sampling 
time being 10 lis. Below: the zoom of the ADC output to show the sign of the first value 
above the background, in this case negative. 



shot at low temperature; in the lower part of the figure the data are zoomed 
around the hit time in order to exhibit the sign of the first swing above the 
background (negative in the shown example); Fig. [6] shows the output in volt of 
the filtering procedure that selects the signal at the first longitudinal resonance. 
The filtered output is shown vs time from the start of the run. 

Fig. [7] shows the Fourier power spectrum for the average of seven signals 
like the one in Fig. [5j after subtracting the noise spectra. The first two lines 
are the first flexural mode and the first longitudinal mode. We have studied 
in detail only the first longitudinal mode, because it is the one of interest for 
gw detectors. However, a quick analysis of the other modes did not show any 
relevant difference regarding the temperature dependence. 

The Fourier component at the desired frequency /o (angular frequency ujq) 
is extracted with a filtering algorithm known as " digital lock-in" . For a hit at 
time corresponding to sample iq, we create the time series: 
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Signal processed with the software iock-in 
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Figure 6: The output of the filtering procedure that selects the signal component at the first 
longitudinal resonance vs time from the start of the run. The shot is the same ofFig.\5\ The 
filtering procedure is most useful for small signals, when the amplitude is comparable to the 
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Figure 7: The Fourier power spectrum for the average of seven signals like that shown in 
Fig. [5[ after subtracting the noise spectrum. The first two spectral lines are the first flexural 
mode and the first longitudinal mode. 
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Cfofeff = ^ V;;cos(a;o — ) 

i=io~N ^'^ 
io+N ^ 

Csignal = ^ ViCOs{uJo--) 
■ ■ Jc 

where Vi is the i — th value of ADC output, /c the samphng frequency and the 
number of samples considered is optimized with N = 400000. Similar quantities 
are computed for the quadrature (sine) component, and from the two we con- 
struct the complex amplitude before (noise) and after (signal) the hit. Taking 
the difference of these two amplitudes produces the desired filtered output. A 
correction is applied to take into account the decay of the signal as function of 
the time. 



4. The RAP measurement in superconducting state 

We present in this section data of the last run of RAP at the Frascati BTF, 
that took place between June 30th and July 2nd 2009, just after the commis- 
sioning of the dilution refrigerator. We took data in the temperature range 
0.344 - 257 K for a total of 164 beam shots on the bar and the energy deposited 
by each shot was in the range 1 ^ ^ 70 mJ. The data at high temperatures 
were in agreement with those of previous run and will not be discussed in this 
section. The 2007 measurements [31], performed at T > 0.54 K, are analyzed 
in this section together with the 2009 data. Some improvements introduced in 
this run, like the "lock-in filter", reduced the noise and increased the sensitivity 
at small energies. 

The data taken at temperatures in the range 0.9 < T < 2 K (i.e. above 
Tc) and T < 2 K, reported in Fig. |8j show a linear correlation between i3, the 
maximum amplitude of the first longitudinal mode, and the released energy W , 
in agreement with the model of section [2] The data of 2009 appear to be of 
better quality than those of 2007, due to the improvements in the beam stability 
and in the analysis procedure. 

Both sets of data for temperature T < 1.6 K of the years 2007 and 2009 are 
shown in Fig.[9j This plot shows the measured B/W (with sign) as function of 
the temperature and of the deposited energy W. The most relevant features of 
this plot are: 

• a constant value B/W for T > Tc, 

• a change of sign oi B/W for T < Tc and 

• a nonlinear dependence of -B on for T <Tc, not predicted by the model 
(Eqn. 
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Figure 8: The measured Amplitude of the first longitudinal mode B^^P vs the energy W 

deposited in the bar in the n state with 0.9 < T < 2 K. The solid line is a linear fit constrained 

to the origin. The top figure represents all data from the 2007 and 2009 runs, while the bottom 

one shows only the 2009 data. The result of the fits are: 

po = (2.33 ± 0.02) X 10-10 [m/J], ^jndf = 51.69 / 55 (2007+2009 plot) 

po = (2.24 ± 0.05) X 10-1° [m/JJ, x^/ndf = 0.4554 / 15(2009 plot) 

The 2009 data have a better due to the improvements in the beam stability and in the data 
analysis. 



A change of sign for superconductive aluminum is to be expected, because the 
effect due to the s n transition can lead to a negative sign due to the com- 
petitive terms in Eqn. ([6| . We recall that the sign of B is inferred by the sign of 
the first value of the ADC over the noise after the beam shot (see Fig. [5]d). The 
sign is positive for an expansion, negative for a contraction. The plot includes 
the measurement of NAUTILUS with cosmic rays described in section [sj This 
value has been obtained from the ratio of the NAUTILUS data at T = 0.14 K 
and the NAUTILUS data at T = 2 K. The NAUTILUS point correspond to a 
value of the energy W ~ 0.5 fiJ. 

Fig. [To] shows the projection of these data in 4 intervals of deposited energy 
W. This figure helps to understand the data behavior, but is important to note 
that some time data don't overlap well, due to the strong dependence on W. 
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B/W vs W and T- data 




Figure 9: Synoptic view of the data for temperature T < 1.6 K, the transition temperature is 
about 0.9 K. The plot shows the measured B/W (with sign) vs temperature T and deposited 
energy W . The most relevant feature of this plot are: a constant value of B/W for T > Tc, 
the change of sign of B/W for T < Tc and the dependence on W of B/W for T < Tc- The 
experimental data are the open circles. The shadowed regions are interpolations of the data. 
The point at the lowest temperature T = 0.14 is obtained from the cosmic ray NAUTILUS 
data. 



The error bars are due to the combination of different sources; 

• the noise in the measurement of the vibration amphtude: ± 1.3 x 10^^'^ m, 

• the uncertainty on the deposited energy (due to a reading error of the 
beam current): ±0.8 mJ, 

• the error in the temperature measurement, that also takes into account 
the local increase in the bar temperature after every shot and a possible 
non uniform profile of the temperature along the bar: ± 0.01 K. 

• an overall systematic error of the order of ± 6%, that accounts for the 
slightly different set-up and analysis procedures adopted in the 2007 and 
2009 runs. 



Fig. 10 shows that the data at the lowest deposited energies have a simple 
behavior: a plateau at very low temperatures, a plateau at higher temperature 
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Figure 10: Summary of the data for temperature T < 1.1 K. The plots show the projections 
of the data of Fig. ^in 4 energy interval. Note that, due to the strong dependence on W , 
data points sometimes do not overlap well. 



(above Tc) and a transition region in between. The first plateau disappears in 
the plots for higher energies. 

As our investigation is aimed at understanding the interactions of cosmic 
rays with a gw detector, we need a model to make prediction of B/W at very 
small value of W: we have used the model described in section [2] adding to it, 
as suggested by the data, a possible saturation of the s ^ n transition effect, 
due the high energy density in the volume crossed by the beam. 

We can estimate this saturation effect starting from the radius of the cylin- 
drical volume that switches from the s to n state around a particle leaving 
195 MeV in the RAP aluminum bar (195 MeV is the mean value of the energy 
loss per each 510 MeV electron of the BTF beam). This critical radius Rc de- 
pends from the energy necessary to activate the transition. This energy can be 
computed in several ways [15j, yielding consistent results. Using Eqn. ([s]) and 
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^ we have: 



where I ^ 2R is the length of the path of the electron inside the bar. 

Therefore, using the value of specific heat of A15056 previously discussed and 
reported in ^24j, we obtain Rc ^ 1 fira ii T — 0.5 K. The total cross surface in- 
terested by this transition for N = 10^ electrons is of the order NttR^ ^ 30 cm^, 
comparable to the beam cross section of the BTF beam, typically ^ 20 cm^. 
This is a crude estimate, but it shows that this effect can produce a non linear 
(with respect to the deposited energy W) response in the RAP data. 

In order to extrapolate the RAP results to values of W and T outside the 
measured range, we used the following four parameters fit to the data for T < Tc'. 

f ^« + (5(T)-«)exp(--^) (9) 

b{T) = + P2T + p3T^ (10) 

Here a ~ 2.25 x 10"^° m J^^ is the constant value of B/W for T > 
obtained from Fig. [s] and b{T) the value of B/W for T < and ^ 0. 
b{T) ^ — 10~^ m J for T = 0.5 K is a function weakly dependent on T and 
accounts for small variations of physical parameters at low temperatures. C/ is 
the integrated specific heat between T and the critical temperature, as defined in 
Eqn. Q and computed from the numerical values of Ref. • Eqn. ([9| derives 
from the consideration that if an electron crosses a region that has already 
undergone the s — n transition, the response is the one of the n stat^ The 
parameter po is the bar volume intercepted by the electron beam. 

In the measurement on a niobium bar |32| the complex pattern of Fig. [9] 
with the non linear behavior in W, was not observed. This is consistent with 
our model if we consider that niobium has a higher Tc and specific heat is much 
larger: C/ of niobium is indeed about two order of magnitude larger than C/ of 
aluminum for comparable value of the integration interval Tc — T. Therefore the 
beam intensity in the niobium measurement was not enough to see saturation 
effects. 

The result of the fit of Eqn. ([9| with these 4 parameters are given in Table 2. 
The data fit gives a /d-o.f. = 368/286 = 1.29, slightly larger than one; this 
suggests that there are effects not taken into account. 



^Consider a beam bunch of high intensity, lasting a few nanoseconds (a time much shorter 
than the relaxation time needed to restore the s state) and involving a volume Vb of the bar: 
the electrons that follow will probe a variable fraction of the volume Vb switched from the s to 
the n state. The volume dV that undergoes the the s n transition for a deposited energy 
dW is 

V VoJ pCi ^ ' 
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parameter 


value 


Po [m^] 
Pi [m J"^] 
P2 [m J-i K-i] 
P3 [m K-2] 


(1.88 ±0.06) X 10-3 
(0.99 ±0.13) X 10-3 
(-1.31 ±0.43) X 10-3 
(-3.0 ±0.37) X 10-3 



Table 2: Parameters for the fit of Eqn. |^ to the RAP data. 



B/W fit residuals 
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Figure 11: Residuals of the fit of Eqn. ^ to the RAP data vs temperature, overlapping on 
all energies. 



Indeed, this simplified transition model cannot include all the effects: for ex- 
ample, the beam profile, the shower development inside the bar, the uncertainty 
of the critical temperature, possible unhomogeneity of the material, etc. are not 
accounted for. It is however remarkable that a very simple expression can de- 
scribe the complicated pattern of Fig. [9j 



Fig. 11 shows the fit residuals (measured data minus the fit predictions) vs 



temperature. The energy dependence seems to be well reproduced with the 
exception of the region just below where there are more scattered points. 

Table 3 summarizes the RAP results of B/W for A15056 obtained with the 
RAP data from room temperature down to 0.14 K. The value and the errors 
for T < Tc are obtained from the fit of Eqn. (|9| for — )• 0. The last value 
was obtained by extrapolating also with respect to the temperature down to 
T — 0.14 K, a temperature relevant for the comparison with NAUTILUS data 
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Temperature 
(K) 


Experimental B/W 
[IQ-io m J-i] 


Predicted B/W 
[10-1" m J-i] 


264 


2.14 ± 0.13 


2.23 


71 


2.18 ± 0.13 


2.23 


4.5 


2.09 ± 0.13 


1.80 


1.5 


2.25 ± 0.13 


1.72 


0.8 


-0.95 ± 2.4 


-7.3 


0.7 


-4.2 ± 2.5 


-7.9 


0.6 


-6.9 ± 2.5 


-7.7 


0.35 


-10.8 ± 2.3 


-9.1 


0.14 


-11.1 ± 1.8 


n/a 



Table 3: Summary of the values of B/W derived from measurements on A15056. ForT > Tc 
the table shows measured data with errors due to the 6% systematic error. Both values and 
errors for T <Tc are obtained by the fit in the low energy limit W 0, taking into account 
the error matrix. The last value at T = 0.14 K was out of the RAP experimental reach, 
and is derived by extrapolating the fit of Eqn. Q also with respect to temperature. The right 
column shows the predictions of the model (see section^^ . 



that was, unfortunately, beyond the reach of our refrigerator. We note that the 
model described in section [2] is quite accurate only for T > Tc- For T < Tc there 
are the discrepancies that could be due either to a failure of the model or to 
uncertainties in the A15056 superconductive parameters. 
Using this table we can estimate: 

X"^°"" =4.9±0.8 (12) 

in good agreement with the value 4.3 ± 1.5 obtained comparing cosmic rays 
detected in NAUTILUS at T = 0.14 K and T 3 K. We assume in the next 
section that this value doesn't change in the range 1.5 ~ 4.5 K. 



5. Cosmic Rays in gw acoustic detectors - interpretation in the light 
of RAP results 

The aim of this paper is to use the RAP results to interpret the cosmic 
ray signals detected in the EXPLORER and NAUTILUS gw antennas. To this 
purpose, we first summarize the most relevant results on this issue [33 | I34 j . 

The ultra-cryogenic acoustic gw detector NAUTILUS [H [3S] is operating 
since 1996 at the INFN Frascati Laboratory, at about 200 meters above sea 
level: it consists of a 3 m, 2300 kg, A15056 alloy bar. We consider here the run 
of 1998, when NAUTILUS was operated at 140 mK. The quantity to be observed 
in this kind of detector (the "gw antenna output") is the vibrational amplitude 
of its first longitudinal mode of oscillation. This is converted by means of an 
electromechanical resonant transducer into an electrical signal which is amplified 
by a dc-SQUID superconducting amplifier. The bar and the resonant transducer 
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form a coupled oscillator system, with two resonant modes, whose frequencies 
were, in that run, /_ = 906.40 Hz and /+ = 921.95 Hz. NAUTILUS is equipped 
with a cosmic ray detection telescope made of seven layers of gas detectors 
(streamer tubes) for a total of 116 counters [3]. 

The gw detector EXPLORER [3il37], similar to NAUTILUS, was located 
in CERN (Geneva-CH) at about 430 meters above the sea level. Scintillator 
counters were installed at EXPLORER in 2002, using scrap equipment recovered 
after the LEP shutdown. The two detectors have a long record of coincidence 
runs [38], also with other detectors [39], to search for gw signals. 

The signal expected in a gw detector like NAUTILUS, as a consequence of 
the interaction of a particle releasing an energy W is [H 1301 SI] , according to 
the model described in section [2] is {W in GeV units): 

£; ^ !|1 X 10^9 5^ [K] (13) 

where the bar oscillation energy E is expressed, as usual in the antenna jargon, 
in Kelvin units (1 K = 1.38 x 10"^'^ J), the numerical constant is the value 
computed using the thermal expansion coefficient and the specific heat of pure 
aluminum at 4 K and (5 is a parameter that describes the difference respect to 
pure aluminum at 4 K. In the previous section we have shown that RAP has 
measured (5„ = 1.16 above the s transition temperature or 5s — 5.7 (= 4.9 x 1.16) 
for superconductive A15056. The vibrational energy E of the first longitudinal 
mode of oscillation is proportional to the square of the Amplitude (in an over- 
simplified model, E = \Muj'^B'^). The constant 7.64 x 10^^ applies if the 
energy is released in the bar center. If the energy is uniformly distributed along 
the bar, as in the case of extensive air showers (EAS), this value is reduced by 
a factor 2. 

Under simplified approximations on the development of the electromagnetic 
shower in the bar, we can derive [H |40l Hi] the relation between the vibrational 
energy detected in the bar and the density A of secondaries in the shower: 

£; = 4.7 X 10^1° A^ 6^ [K] (14) 



The plot of vibrational energy E^xp vs particle density A is shown in Fig. 12 
In this figure we show the events detected by NAUTILUS, both at 140 mK and 
at 2.6 K, as well as by EXPLORER at a temperature of about 3 K. We clearly 
see a difference of more than one order of magnitude between the measurements 
taken with aluminum in the s state and those in n conduction state. From this 
plot we can estimate a value of 6s — 5.0 ± 1.8 in good agreement with the value 
derived by the RAP experiment: 6s = 5.7 ± 0.9. 

We now briefly discuss the event rate of cosmic rays in gw detectors. The 
cosmic ray event rate in NAUTILUS and EXPLORER has been evaluated in 
the past considering three different event categories: i) pure electromagnetic 
showers, responsible for most of the high energy events detected in the bar 
detectors; ii) showers produced by muons and iii) showers produced by hadrons 



in the bar. We use Eqn. (13) with the correction 6n = 1.16 for the response of 
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Figure 12: Averages of signals with energy E^xp < 0.1 K, grouping data in ranges of particle 
density A. Filled circles NAUTILUS at T = 0.14 K, open circles NAUTILUS atT = 3 K, 
filled squares EXPLORER at T = 3 K. The data gathered at T = 0. 14 K are roughly one 
order of magnitude larger than those collected at T = 3 K, From Ref. \34l - 



an A15056 bar in the n state. The rate of electromagnetic air showers (EAS) 
is computed starting from the empirical relation due to G. Cocconi [12] ■ The 
event rate due to muon and hadron interactions inside the bar was computed 
using the GEANT package [43], to simulate the antenna and the CORSIKA 
Monte Carlo [H], as input to GEANT, to simulate the effect of the hadrons 
produced by the cosmic ray interactions in the atmosphere, assuming a cosmic 
ray "light" composition. The Monte Carlo simulation represents 1 year of data 
taking. 

The results are shown in Table g] The energy in the first longitudinal mode 
E (first column) is proportional to the square of the absorbed energy W . 
The rate of the the events scales as W~^ '^. This is because the cosmic ray 
integral spectrum is well described by a power law E'^^ with /3 ^ 1.7 for cosmic 
ray primaries up to the so called "knee" at Ecr — 10^^ eV and /? ^ 2 at higher 
energies. 

In Ref. [40], very large NAUTILUS signals at a rate much greater than 
expected were reported. In the light of the analysis reported above, it is now 
clear that the value 5s = 5.7 must be used in Eqn. (14) to compute the expected 
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Vibrational E 


Deposited W 


Total 


[K] 


[GeV] 


[events/day] 


> 10-^ 


> 44.5 


107 


> 10-* 


> 141 


14.5 


> 10-3 


> 445 


1.6 


> 10-2 


> 1410 


0.19 


> 10-1 


> 4450 


0.03 



Table 4: Estimated rate of antenna excitations due to cosmic rays in NAUTILUS vs the 
vibrational energy of the longitudinal fundamental mode that such events can produce. The 
value at E = 0.1 K is obtained extrapolating from the lower energy values. The values in the 
second column are the values of cosmic ray energy that the bar needs to absorb in order to 
have an excitation energy E. Vibrational and Deposited energy are correlated by Eqn. 
under the assumption of uniformly distributed energy. 



response and the expected rate. The NAUTILUS 1998 data event rate per 
day after the unfolding of the background, with the procedure described in 
[34] is shown in Fig. 13 The continues line is the prediction of Table |4] with 
6s — 5.7. We find now a good agreement between measurements and predictions; 
previously hypothesized exotic explanations, based on anomalous component of 
cosmic rays or anomalous interactions of cosmic rays with a superconductive 
bar can now be excluded. 

Cosmic ray showers are a very important tool to verify the sensitivity of gw 
bar detectors to signals distributed along the bar, signals similar to gw even if the 
excitation mechanism is different. Moreover the cumulative analysis of Fig. [12] 
and the agreement with the RAP measurements show that for this kind of search 
the sensitivity of bar detectors can be studied down to 10 /xK corresponding to 
B < 10~i^ m. This results is quite important for other "cumulative" analysis 
like the search of signals from gamma ray bursts. 

Finally, we remark that acoustic gw detectors have no limitations due to 
saturation effects in detecting large signals. Indeed the largest event detected 
up to now has a vibrational energy in the first longitudinal mode E ~ 670 K 
corresponding to ^ 360 TeV in the bar. The event occurred in EXPLORER on 
Nov 10 2006 9:40 UT. 



6. Conclusions 

We have shown that the Thermo- Acoustic Model reasonably well describes 
the response of a bar to the passage of ionizing particles. At high temperatures, 
in a normal conduction state, the prediction, only based on the knowledge of 
the Griineisen parameter, is in a good agreement with the data. In a super- 
conductive material, the transitions between the s and n state complicate this 
picture. Due to the poor knowledge of the low temperatures parameter and to 
the approximations of the model, a direct measurement with a particle beam 
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Figure 13: NAUTILUS 1998, atT = 0.14 K. The integral distribution of the event rate after 
the background unfolding, compared with the expected distribution (continuous line). The 
prediction is computed using the data of Tableland using the value Ss = 5.7 measured by 
RAP. The good agreement suggests the absence of anomalous components of cosmic rays or 
anomalous interactions of cosmic rays with a superconductive bar. Modified from Ref. ^3^. 



was needed to directly measure the response to ionizing particles. The RAP 
experiment addressed this issue. At high energy densities of the impinging par- 
ticles, we have detected and studied with RAP non linear effects that complicate 
the data analysis. 

We have shown that the unexpected large events detected in 1998 with NAU- 
TILUS at T = 0.14 K were due to its superconductive state. Using the RAP 
measurements we have regauged the rate of cosmic rays detected by the NAU- 
TILUS and EXPLORER antennas, and shown that they are in agreement with 
the predictions. 

Currently the background due to cosmic rays in acoustic bar detectors is 
negligible. This is because the typical sensitivity \s E = 1 mK; the standard 
event selection requires a threshold of about 25 mK. With such a threshold we 
have a few events per week due to cosmic rays. These events are however very 
useful as a tool to continuously monitor and calibrate the acoustic gw detectors. 
Moreover, in the standard NAUTILUS data analysis, these events are currently 
vetoed and removed by the official gw event list. 

Cosmic rays will represent an important source of background in future 
higher sensitivity [7] , possibly superconducting, detectors, and this noise should 
be taken into account in both acoustic [151 H51 WJ\ and interferometric [IS] 
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detectors. To remove this background, moving to an underground site could be 
necessary. 
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